function logpsilabor = loglaborprod(s,x,p)

% note: if x and s are vectors, the function returns
% [psi(s(1),x(1)), psi(s(2),x(2)),...,psi(s(end),x(end))]
% x and s must have the same length

logpsilabor = p.m*s + p.n*s.^2 + p.a*x.*s + p.A;

end